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I. INTRODUCTION 



Understanding properties of hadrons directly from QCD represents a considerable on- 
going challenge. A feature of strong interaction dynamics is the confinement of quarks and 
gluons into color-neutral hadronic states, and their spectrum is a fundamental observable of 
the theory. Calculation of hadronic masses will be a milepost in strong interaction physics. 
For baryons containing b or c quarks, the underlying dynamics encompasses the interplay of 
heavy-quark physics and light-quark physics, which are opposite limits of QCD. 

Lattice QCD provides a first principles tool for numerical determination of QCD ob- 
servables, for example, the mass spectrum of baryons containing a heavy quark. Cal- 
culations of singly heavy baryon masses have been pursued in quenched QCD (QQCD) 
QCD 0, 0, H 0, IE El using heavy quark effective actions. Up until now 0], partially 
quenched QCD (PQQCD) calculations in the baryonic sector have been very limited. The 
necessity for quenched and partially quenched variations approximations in lattice QCD 
arises from computational difficulties inherent to calculating the fermionic determinant. In 
QQCD, the time-costly determinant is replaced by a constant, whereas in PQQCD the de- 
terminant is calculated by using larger quark masses for the quarks not coupled to external 
sources. In PQQCD, sea quark contributions are thus retained; while in QQCD, the sea 
quarks have been unphysically integrated out. Due to restrictive computational time, the 
valence quark masses, too, cannot be set to their physical values and one must extrapolate 
down from the values used on the lattice. 

To perform rigorous extrapolations in the quark mass, one must understand how QCD 
observables behave as the quark masses vary. For low-energy QCD, there exists a model- 
independent effective theory, chiral perturbation theory (%PT), that is written in terms 
of the pseudo-Goldstone bosons appearing from chiral symmetry breaking 0, While 
these bosons are not massless, they remain light compared to the scale of chiral symmetry 
breaking and hence dominate low-energy hadronic interactions. Because the light quark 
masses explicitly break chiral symmetry, they appear as parameters of and enable 
the determination of the quark mass dependence of QCD observables by matching onto 
the effective theory. Without external input (e.g. experimental data or lattice QCD data), 
the effective theory generally lacks complete predictive power as symmetries constrain the 
types of operators in the Lagrangian, not the values of their coefficients. While providing 
a guide to extrapolate lattice QCD data, %PT also benefits from the determination of its 
low-energy constants (LECs) from these data. Furthermore heavy quark symmetry can 
be combined with chiral symmetry to describe hadrons formed from both heavy and light 
quarks U , 12 • 



For QQCD lattice simulations, such as those for singly heavy baryons, quenched chi- 
ral perturbation theory (QyPT) has been developed for the extrapolation in valence quark 
masses [l3L FLU EH EH fl7| . There is, however, no general connection of QQCD observables to 
QCD because QQCD lacks an axial anomaly. As a manifestation of this fact, Qx?T contains 
operators that do not have analogues in %PT. Furthermore, the LECs of the operators that 
have analogues are numerically unrelated to their counterparts in %PT. The quenched ap- 
proximation can be remedied by performing partially quenched lattice simulations. For such 
simulations, partially quenched chiral perturbation theory PQ^PT has been constructed to 
perform the extrapolation in both sea and valence quark masses 

Unlike QQCD, PQQCD retains an axial anomaly and the flavor singlet field can be inte- 
grated out. Consequently the LECs of xPT appear in PQxPT, and hence there is a means 
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to determine %PT parameters using PQQCD lattice simulations. Much work has resulted 
from this possibility, e.g., in the baryonic sector see 0, H H H m, H H H3 . 



We anticipate the inclusion of effects from dynamical light quarks into the singly heavy 
baryon sector, and in this work we determine the masses of baryons containing a heavy quark 
at next-to leading order (NLO) in partially quenched heavy hadron xPT. Throughout we 
work in the heavy-quark limit. Higher-order corrections can be considered if one wishes to 
make contact with simulations that go beyond the leading-order heavy quark effective action. 
We determine the masses for three-flavor and two-flavor partially quenched simulations and 
use the isospin limit in the former. Such expressions allow one to perform the required 
extrapolation in light-quark masses. Additionally the LECs appearing at NLO in heavy 
hadron xPT can be determined. We spell out explicitly the LECs that appear in the mass 
splitting of the S£ = and S£ = 1 baryons from chiral effects to this order, and investigate 
the chiral contributions to the mass relation among the sextet baryons. 

The organization of the paper is as follows. First in SeclHl we include baryons containing 
a single heavy quark into PQ;\PT. We write down the terms of the SU(Q\3) heavy hadron 
chiral Lagrangian that are relevant for the determination of baryon masses at NLO. The 
masses of the se = baryons are calculated in Sec. IHIl while those of the se = 1 baryons 
are calculated in Sec. IIV1 The corresponding results for S77(4|2) with non-degenerate light 
quarks are presented in Appendix EI For reference, we give the masses in SU(3) and SU(2) 
%PT in Appendix |B] The summary, Sec. EJ highlights the goal of determining the mass 
splittings between baryons containing a single heavy quark. 



II. PQxPT FOR HEAVY HADRONS 

In this section, we formulate PQxPT for heavy hadrons in SU(6 | 3). We begin by review- 
ing PQQCD. Next we review the pseudo-Goldstone mesons of PQ^PT, and finally include 
the heavy hadrons into this partially quenched theory. 



A. PQQCD 

In PQQCD, the Lagrangian in the light-quark sector is 

9 

j,k=i 

This differs from the usual three light-flavor Lagrangian of QCD by the inclusion of six 
extra quarks; three bosonic ghost quarks, (u, d, s), and three fermionic sea quarks, (j,l,r), 
in addition to the light physical quarks (u,d,s). The nine quark fields transform in the 
fundamental representation of the graded SU(6\3) group 3^, 3- They have been placed in 
the nine- component vector 



q = {u,d,s,j,l,r,u,d,s) T . (2) 
These quark fields obey graded equal-time commutation relations 

<?f (x)gf(y) - (-l)^gf (y)tf(x) = rt^ 3 (x - y), (3) 



3 



where a, (3 and i,j are spin and flavor indices, respectively. The vanishing graded equal-time 
commutation relations are written analogously. Above we have made use of grading factors 

f 1 for k = 1,2,3,4,5,6 
Vk ~ \ for k = 7, 8, 9 ' W 

that take into account the different statistics of the PQQCD quark fields. When m u = m<f, 
the quark mass matrix of SU(6\3) is given by 

m q = diag(m u , m u , m s , rrij, rrij, m r , m u , m u , m s ). (5) 

Notice that each valence quark mass is identically equal to its corresponding ghost quark's 
mass. Thus there is an exact cancellation in the path integral between their respective 
determinants. The only contributions to observables from disconnected quark loops come 
from the sea sector of PQQCD. Thus one can keep valence and sea contributions separate 
and additionally vary the masses of the valence and sea sectors independently. In the limit 
rrij — > m u and m r — > m s , QCD is recovered. 

B. Mesons of PQxPT 

In the limit that the light quarks become massless, the theory corresponding to the La- 
grangian in Eq. (JTJ) has a graded symmetry of the form SU(6\3)l (g) SU(6\3)ji <8> U(l)y, 
which is assumed to be spontaneously broken down to a remnant SU(6\3)y ® U(l)v sym- 
metry in analogy with QCD. One can build an effective low-energy theory of PQQCD by 
perturbing about the physical vacuum state. This theory is PQxPT, and the emerging 
pseudo-Goldstone mesons of this theory have dynamics described at leading order in the 
chiral expansion by the Lagrangian 

f 2 

£ = ±-str (d^d^) + A str (m ? S f + mJS) + a^%d^% - m^, (6) 
where the field 

S = exp(^)^, (7) 
and the meson fields appear in the SU(6\3) matrix 

The quantities a$ and mo are non- vanishing in the chiral limit. The M and M matrices 
contain bosonic mesons (with quantum numbers of qq pairs and qq pairs, respectively), 
while the x an d matrices contain fermionic mesons (with quantum numbers of qq pairs 
and qq pairs, respectively). The upper 3x3 block of the matrix M contains the familiar 
pions, kaons, and eta, while the remaining components consist of mesons formed from one 
or two sea quarks, see e.g. j2^|. The operation str() in Eq. (jUJ) is a supertrace over flavor 
indices. Expanding the Lagrangian Eq. (jUJ) to leading order, one finds that mesons with 
quark content qq 1 are canonically normalized when their masses are given by 

rn 2 qq , = ^(m q + m ql ). (9) 
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The flavor singlet field that appears above is defined to be $ = str($)/ v / 6. As PQQCD 
has a strong axial anomaly U(1)a, the mass of the singlet field mo can be taken to be on 
the order of the chiral symmetry breaking scale, and thus $o in Eq. © is integrated out 
of the theory. The resulting rj two-point correlation functions, however, deviate from their 
familiar form in %PT. For a,b = u, d, s,j, I, r, u, d, s, the leading-order r) a r]b propagator with 
2 + 1 sea-quarks is given by 



e a t b (q 2 - mf ) (q 



mz 



q 2 — m 2 a +ie 3 (g 2 — m 2 + ie) {q 2 — m 2 b + ie) (q 2 



m x + ie > 



where 



(a 



-1 



(10) 



11^ 



The mass m Vx is the mass of a quark-basis meson composed of x and x quarks, while the mass 
mx is defined as m\ = ~ (m 2 ^ + 2m 2 r ) . The flavor neutral propagator can be conveniently 
rewritten as 

GriaVb = e a$abPa + ^a^b^ab (P a , Pb, Px) i (12) 

where 



Pa 



mi. +ie' 



m^ + ie ' 



9 2 i • ' 

q z — m x + ie 
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C. Baryons containing heavy quarks in PQyPT 



Heavy quark effective theory (HQET) and PQxPT can be combined to describe the 
interactions of heavy hadrons with the pseudo-Goldstone modes. The hadrons we treat 
below are baryons composed of one heavy quark Q, where Q = c or b, and two light quarks. 
Let mq denote the mass of the heavy quark. In the limit mq — > oo, these baryons (and 
indeed all heavy hadrons) can be classified by the spin of their light degrees of freedom, 
S£. This classification is possible because the heavy quark's spin decouples from the system 
in the mq — > oo limit. The inclusion of baryons containing a heavy quark into yPT was 
carried out in 0, |3^, |36| and their masses were calculated to one-loop order in J37| . The 
quenched chiral theory for baryons with a heavy quark was written down in [38( and the 
partially quenched theory was investigated in [27] for SU(A\2). Here we include the baryons 
in S77(6|3) PQxPT. 

Consider first the baryons with si = 0. These have spin | because there is only one 
way to combine the spin of the heavy quark with that of the light degrees of freedom. To 
include these baryons into PQy PT, we use the method of interpolating fields to classify their 



representations of SU(6\3) jl7J, |2J, |25|, |27J, |39[. The si = baryons are described by the 
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field 1 



'ij 



Q 



7,c 



a, a 8,b , 8,b a,a\ / \ 

Qi Qj + Qi Qj ) £abc{C7 5 )a(3, 



(14) 



where i and j are flavor indices. The flavor tensor T has the property 



T 

1 i.i 



(15) 



and thus forms a 39-dimensional representation of SU{6\3). It is convenient to classify these 
states under the quark sectors (valence, sea or ghost) acted upon and to this end we use the 
super-algebra terminology of BE}. Under SU(3) yal ® 5C/(3) sea ® S77(3) ghost , the ground 
floor, level A transforms as a (3, 1, 1) and contains the familiar se = baryon tensor of 
QCD Ty, explicitly %j = T^, when the indices are restricted to the range 1 — 3. In our 
normalization convention, we have 



/ 



Ta 



1 

71 






A n T 5 \ 


-Aq 


o ^ 




/ 



(16) 



where the superscript labels the 3-projection of isospin. The first floor of level A contains 
nine baryons that transform as a (3, 1, 3), while the ground floor of level B contains nine 
baryons that transform as a (3, 3, 1). The remaining floors and levels are simple to construct 
but are not necessary for our calculation. 

Next we consider the S£ = 1 baryons. Such baryons come in two varieties, spin | and |, 
because there are two ways to combine se with the heavy quark's spin. In the uiq — > oo 
limit, these two multiplets are degenerate and can be described by one field S^. When 
establishing the representations, we are of course ignoring this two-fold degeneracy (which 
is only lifted away from the heavy quark limit). The interpolating field for these baryons is 



Q 



7,c 



a, a 8,b 

Qi Qj 



8,b a, a 

Qi Qj 



(17) 



where the flavor tensor satisfies 



CM 

6 ij 



' "if 



(18) 

and hence makes up a 42-dimensional representation of S77(6|3). Under 577(3)^ <g> 
SU(3) sea _ <S> 5'f^(3)ghost, the ground floor, level A transforms as a (6, 1, 1) and contains the 
familiar si = 1 baryon tensor of QCD S^, explicitly = S^, when the indices are restricted 
to the range 1 — 3. Here the QCD flavor tensor is given as 



CM 



7 ^) 7 5 By + Bg, 



(19) 



with 



B 



ij 



Y+ 1 
Q 

l ^0 
Q 



1 vO 

s^ 1 -4 



n 



(20) 



IJ 



1 Since we shall work only to leading order in 1/mg, we omit the label Q from all baryon tensors. 
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and similarly for the Above v M is the velocity vector of the heavy hadron and we have 
suppressed velocity labels on all heavy hadron fields throughout. Again the superscript on 
these states labels the 3-projection of isospin. The first floor of level A contains nine baryons 
that transform as a (3, 1, 3), while the ground floor of level B contains nine baryons that 
transform as a (3, 3, 1). The remaining floors and levels are simple to construct but are not 
necessary for our calculation. 

The free Lagrangian for the T and fields is given by 

£ = -i (S»v ■ DS^) + A (S"5 M ) + Xt (S^M+S^) + A 2 (S^) stiM + 

+i (Tv ■ DT) + A 3 (TM + T) + A 4 (TT) strM + . (21) 

Above we have employed Q-notation for the flavor contractions that are invariant under 
chiral transformations. The relevant contractions can be found in [2?}. The mass of the T 
field has been absorbed into the static phase of the heavy hadron fields. Thus the leading- 
order mass splitting A appears as the mass of the iS M . This splitting remains finite as 
rriQ — > oo and cannot be removed by field redefinitions due to the interaction of the T and 
iS M fields. In our power counting, we treat A ~ m n . The Lagrangian contains the chiral- 
covariant derivative D^, the action of which is identical on T and fields and has the 
form [27^ 

{D^T\ 3 =d^ + V^ 3 + {-r^ + ^V^% r . (22) 
The vector and axial-vector meson fields are defined by 

V = \ + $r£) , A»= l - (£c^t - £t«n) , (23) 

respectively. The mass operator appears as 

M+ = \ {im q i + ^m q £}) . (24) 

To calculate baryon masses to NLO, the above Lagrangian is not sufficient at tree level. 
This is because the small parameter A transforms as a singlet under chiral rotations. Con- 
sequently arbitrary polynomial functions of A/A x must multiply terms in the Lagrangian, 
for example one such term is 

We will not write down all operators that lead to these effects nor will we keep track of 
similar such terms that arise from loop integrals. To account for such terms to the order we 
are working, we must treat the LECs as linear functions of A. Specifically 

A^A 3 (A). Af+A <^. (25) 

Determination of the LECs for the four omitted A/A x operators requires the ability to vary 
A, and for this reason we keep the effect of these operators only implicitly. 2 



One must also consider operators that change the value of the mass splitting in the chiral limit. Let us 
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FIG. 1: One-loop graphs which give contributions to the masses of the sg = baryons. The single, 
double and dashed lines correspond to sg = baryons, sg = 1 baryons and mesons, respectively. 
The filled squares denote the axial coupling given in Eq. (|26|). while the cross denotes the hairpin 
interaction. 



The Lagrangian that describes the interactions of the T and S 11 fields with the pseudo- 
Goldstone modes is given at leading order by the Lagrangian 

£ = ig 2 (SWAPS*) W + V2g 3 [(TA^) + (5"i,T)] , (26) 

and gives rise to the loop corrections to the baryon masses at NLO. The LECs appearing 
above, the \j, g 2 , and g%, all have the same numerical values as those used in 577(3) heavy 
hadron x?T, see Appendix IB1 



III. MASSES OF THE s e = BARYONS 



The mass of the i = baryon, Tj, can be written in the form 

M Ti = M (/x) - Mg> 0*) - Mg /2) (a*) - ... . (27) 

Here, M (/x) is the renormalized mass of the T baryons in the chiral limit. This quantity 
is independent of m q and also of the Tj. My. is the contribution to the i th si = baryon 

mass of order rriq, and /x is the renormalization scale. In this calculation we use dimensional 
regularization with a modified minimal subtraction scheme, where we have subtracted terms 
proportional to 

- -j E + l + log47r. (28) 



denote A^ ) as the mass of the iS M field in Eq. 121|) . To the order we are working, we need the following 
operators 



£ = AW 



S^) + A«^ (S^S,) + (S»S, 



An overall factor of A^ ^ appears to render the new LECs A^ and A^ 2 ^ dimensionless. The mass splitting 
A that results from these contributions is 



A = A<°> 



A x V Ax / 



Additionally the loop integrals contribute to A up to 0(A^ ). Instead of treating A as a function of 
A (o) 

we combine the new LECs and contributions to A from loop integrals into the parameter A that 
appears in Eq. J2U- Treating A as the mass of the field instead of A^ ^ results in differences that are 
of higher order than we need consider. 
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TABLE I: The tree-level coefficients in xPT and PQxPT. The SU(3) and SU(Q\3) coefficients tut 
are listed for the si = baryon states T. 



T 


rriT 


Aq 


m u 


E Q 


\{m u + m s ) 



TABLE II: The coefficients and A^, in SU{6\3) PQxPT. Coefficients are listed for the s £ = 
baryon states T, and for A^ are grouped into contributions from loop mesons with mass m^, while 
for A T ,,i are grouped into contributions from pairs of quark-basis r\ q mesons. 
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To calculate the masses of the = baryons in PQxPT, we need the tree-level contri- 
butions from the operators in Eq. (|21|). These give rise to the leading-order contributions 
to the mass, namely 

jl4 1} = X 3 m T + A 4 str(m 9 ), (29) 

where the coefficients ttit are listed for individual T-states in Table [I] 

There are also loop diagrams which are generated from the vertices present in Eq. (|2fi|) . 
The relevant loop contributions are depicted in Fig. Q and include the hairpin interaction 
from the flavor-neutral field. These two one-loop graphs contribute to the T masses at NLO, 
explicitly we find 



M, 



(3/2) 



9l 



4vr 2 / 2 



(30) 



The non-analytic functions T appearing in the above equation are defined by 

5 — -\A 2 — m 2 + is 



T{m, 5, /i) 



[m 



8 2 



V<5 2 — m 2 log 



5 + V<5 2 — rn 2 + is 



— -5m 2 log 



m 



and 



(31) 



T =Hw[J 7 (m$,5,tx),J r (m 4>/ ,5,ij),J r (mx,5,ijL)], (32) 

for the hairpin contributions. The first sum in Eq. (J5U|) is over all loop mesons of mass 
m^. The second term sums over all pairs of flavor-neutral states in the quark basis, e.g., 
above (fxp' runs over r] u r] u , r] u 7] s , and rj s 7] s ; and all other entries are zero. In this way there 
is no double counting. The coefficients AT, and AT,, appear in Table ILD and depend on the 
particular baryon state T. 
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TABLE III: The tree-level coefficients in X PT and PQxPT. The SU(3) and 5J7(6|3) coefficients 
m B are listed for the sg = 1 baryon states B. 
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m B 




m u 


e 'q 


\{m u + m s ) 




m s 





FIG. 2: One- loop graphs which give contributions to the masses of the S£ = 1 baryons. The single, 
double and dashed lines correspond to = baryons, s% = 1 baryons and mesons, respectively. 
The filled squares denote the relevant axial couplings given in Eq. (|26|). while the crosses denote 
hairpin interactions. 



IV. MASSES OF THE s e = 1 BARYONS 

The mass of the i th se = 1 baryon, Bi, can be written in the form 

M Bi = M (/x) + A(/i) + M« 0*) + M| /2) (/x) + . . . . (33) 

Here, Mo 00 is the renormalized mass of the T baryons in the chiral limit and A 00 is 
the renormalized splitting between the T and B baryons in the chiral limit. Both of these 
quantities are independent of m q and also of the B^. M B n ^ is the contribution to the i th 

S£ = 1 baryon mass of order and fi is the renormalization scale. 

To calculate the masses of the sg = 1 baryons in PQ^PT, we need the tree-level contri- 
butions from the operators in Eq. (j21|) . These give rise to the leading-order contributions 
to the masses, namely 

M ( B 1] = \ imB + A 2 str(m g ), (34) 

where the coefficients m B are listed for individual 5-states in Table ILTTl 

There are also loop diagrams that contribute to the masses. These are generated from 
the vertices present in Eq. ()26|). The relevant loop contributions are depicted in Fig. El 
and include hairpin interactions from the flavor-neutral field. These four one-loop graphs 
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TABLE IV: The coefficients Aj and A^, in SJ7(6|3) PQxPT. Coefficients are listed for the sn = 1 
baryon states B, and for A^ are grouped into contributions from loop mesons with mass m^, while 
for A?,, are grouped into contributions from pairs of quark-basis rj q mesons. 
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TABLE V: The coefficients Bf and B^, in St7(6|3) PQxPT. Coefficients are listed for the a t = 1 
baryon states B, and for B^ are grouped into contributions from loop mesons with mass m^, while 
for B?,, are grouped into contributions from pairs of quark-basis r] q mesons. 
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contribute to the i? masses at NLO, explicitly we find 



M 



(3/2) 



127T 2 / 2 

\2tiP 



Yl A %^ m ^ - A, /i) + ^ A %^{^ "V '» ~ A ' /") 



(35) 



The non-analytic functions JF appear in Eqs. (J31j) and (|32j) . The function .M arising from 
hairpin contributions is given by 



(36) 



The coefficients A^ , and A^,, appear in Table HVl and depend on the particular baryon state 
B. The values of B? and B?,, are similarly listed in Table IVl 



V. SUMMARY 



Above we have determined the masses of both the Sf = and sg, = 1 baryons at one-loop 
order in the chiral expansion. We performed these partially quenched calculations in the 
isospin limit of SU (6\3) and away from the isospin limit in S77(4|2). Understanding the 
light-quark mass dependence of these baryons' masses in and PQ%PT is necessary to 
extrapolate QCD and PQQCD lattice data from the valence and sea quark masses used in 
the numerical simulations down to their physical values. Because the LECs of are all 
contained in PQ^PT, fitting PQQCD lattice data to PQxPT expressions can yield rigorous 
predictions for QCD. 



11 



The lattice data, moreover, can be used to determine the LECs appearing in the chiral 
theories. At 0(m^ 2 ) in the chiral expansion (and in the heavy quark limit), there are 
only a few LECs which enter into the calculation and lattice data can be used to fit these 
universal parameters. Of particular interest for study on the lattice are the properties of the 
Qq baryon formed from two strange quarks and a heavy quark. This state does not decay 
strongly and current lattice technology can reach realistic values for the strange quark mass. 
As with the Q~ 0], the Qq is thus an ideal testing ground for obtaining and extrapolating 
lattice data with dynamical quarks. The expressions derived above allow one to analyze 
the light-quark mass dependence of the Qq mass in partially quenched simulations. Further 
work is needed for other processes. 

An example of the utility of ;\PT is to predict the widths for strong decays of heavy 
baryons from the lattice, e.g. for Eg — > Aqtt. Treating the resonant Eg state directly on 
the lattice is problematic for real-valued Euclidean space simulations. Instead, one can work 
with pion masses m n > 200 MeV for which the Eg is stable, and where one still trusts 
the chiral expansion. 3 This enables the lattice simulations to be performed and one then 
uses the %PT or PQ^PT expressions derived above to fit the LECs from the lattice data. 
Consequently these LECs lead to predictions for the physical Eg resonances. For example, 
their decay widths can be found by knowing the pion mass and the constants /, gs, and A. 
The expression for the Eg width Sm(Ms Q ) in terms of these parameters is 



,2 



3m(M s J = -^(A 2 -m^/ 2 . (37) 

As a further application, the mass splitting between the sg = and sg — 1 baryons due to 
chiral effects can be studied to isolate the various LECs. For the real part of the difference 
between the average Eg and Ag masses in SU (2) flavor, we have 



AM = 3?e(M SQ - M K 

A + ( ^Ai + A 2 + A 4 ) tr(m 9 ) + ^p m - + A ' ^ (35 

While the chiral contributions to individual baryon masses are relatively small (naively they 
are ~ 5%), the mass splitting between the si = and s# = 1 baryons is an order of magnitude 
more sensitive to quark mass variation. 4 To see this variation, we plot the mass splitting 
AM as a function of the pion mass in Fig. |H1 In order to plot this splitting, we have assumed 
the Xj LECs are of order one at the scale /i = 1 GeV, and hence their contribution to the 
splitting can be neglected. We have used the values g 2 = —1.9 and g^ = 1.5 as estimates of 
the couplings. These values come from large N c considerations (ill . 42 |. The parameter A 
is chosen to be 200 MeV, which is near the phenomenological value. 

Considering the sextet baryons in the isospin limit of SU(3), the quark mass dependence 
from chiral loops can be isolated for the linear combination of the masses 5M, given by 

5M = M Eq + M nQ - 2Mn' . (39) 



3 This is a conservative estimate because finite volume effects and heavy quark effects will modify the £q 
decay threshold. 

4 This is also true in PQ^PT, where there is more freedom to vary masses and determine parameters. 
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FIG. 3: Pion mass dependence of the splitting AM between the si = and S£ = 1 baryons in two 
flavor xPT. 

0.20 1 , 1 , 1 , 1 , 1 
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FIG. 4: Mass dependence of the function 5M involving the sextet baryon masses in the isospin 
limit of xPT. The strange quark mass is fixed using the experimental values for the pion and kaon 
masses. The up quark mass varies and the resulting variation of 5M is plotted as a function of the 
pion mass. 

This is because at LO in the chiral expansion SM = 0, and deviations from zero occur at 
one- loop order. 5 Thus this combination of the masses is only sensitive to the couplings #2 
and g$. Using our expressions from Appendix El we find 

6M = m^p S a * si m l + al H™>h - A > y)] > ( 4 °) 

where the coefficient is given by a,,- = 1/4, ax = — 1, and a v = 3/4. Notice there is 
no /i-dependence in SM because of the flavor structure of the coefficient a^, and the Gell- 
Mann-Oaks-Renner relation. Using the same estimates of g% and A as above, we can 
plot SM as a function of the quark masses. This is done in Fig. 0] for the real part of SM. 



5 Additionally in SU(6\3) PQyPT, the tree-level contributions to 5M cancel leaving the one-loop results 
as the leading non-vanishing contribution. For simplicity we do not consider the extra freedom available 
to extract the couplings by varying the sea quark masses. 
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The imaginary part is an order of magnitude smaller. For simplicity, we hold the strange 
quark mass fixed, and vary the up quark mass. Results are plotted as a function of the pion 
mass. Determination of the quark mass dependence of SM provides a clean way to access 
the squares of the couplings g 2 and g 3 from lattice data or partially quenched lattice data. 

Soon lattice QCD calculations will enable an understanding of the properties of hadrons 
in terms of quark and gluon degrees of freedom. A crucial first step will be the determination 
of their masses from first principles. This requires moving beyond the quenched approxi- 
mation and including effects from dynamical light quarks. Current and foreseeable lattice 
calculations will rely upon chiral extrapolations in the light-quark sector. Our expressions 
for the masses of singly heavy baryons in and PQxPT enable such extrapolations. One 
can go beyond the results here and include higher-order terms in HQET or try to address 
lattice artifacts. Such a study would improve our ability to extrapolate lattice QCD data 
on two fronts: the light-quark regime and the heavy-quark regime. 
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APPENDIX A: HEAVY HADRON MASSES IN 5J7(4|2) 

The set up of heavy hadron chiral perturbation theory in the baryon sector for SU(4\2) is 
quite similar to that of S77(6]3) in Sec. |H] The partially quenched heavy hadron Lagrangian 
has been written down in [27]. In this appendix, we briefly review the setup of S77(4|2) 
PQQCD and PQ%PT. We then present the calculation of heavy hadron masses for this 
graded flavor group. 



1. PQxPT for SU(4\2) 

In SU(4\2) PQQCD, the Lagrangian in the light-quark sector is 

6 

£= ^-(ip-m^ft, ( Al ) 

j,k=i 

and differs from the usual two light-flavor Lagrangian of QCD by the inclusion of four 
extra quarks; two bosonic ghost quarks, (u,d), and two fermionic sea quarks, (J, I), in 
addition to the light physical quarks (u, d). The six quark fields transform in the fundamental 
representation of the graded S77(4|2) group. They appear in the six-component vector 

q = (u,d,j,l,u,d) T . (A2) 

The quark fields obey the same graded equal-time commutation relations as in Eq. © but 
with the grading factor rjk now defined by 
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The mass matrix with non-degenerate quarks is given by 



9 



diag(m u ,m d ,mj,rni,m u ,m d ), (A4) 



where the ghost quarks are still degenerate with their valence partners. In the limit rrij —>■ m u 
and mi —>■ m d , QCD with two light flavors is recovered. 

In the limit of massless light quarks, the Lagrangian in Eq. (jAljl has a graded symmetry 
S77(4|2) L <g> SU(4\2) R ®U(l)v, which is assumed to be spontaneously broken to SU(4\2) V <8 
U(l)v in analogy with QCD. The Lagrangian describing the pseudo-Goldstone mesons of 
this theory is identical in form to Eq. (JHJ). The only difference being that the meson fields 
appearing in Eq. (jHJ) are given by $, which is now an SU(4\2) matrix. The upper 2x2 
block of the matrix M contains the familiar pions, see e.g. |25j. The flavor singlet field that 
appears m S77(4|2) is defined to be $ = str($)/ v / 2. As before the mass of the singlet 
field mo can be taken to be on the order of the chiral symmetry breaking scale. In this 
limit, the rj two-point correlation functions deviate from their familiar form in x?T. For 
a,b — u, d,j, I, u, d, the leading-order r\ a r\ b propagator is given by 

G ie a 5 ab i e a e b (q 2 - m^) (g 2 - m|) 

y ™„ q2 _ ^ +ie 2 [f - m 2 Va + ie) (q 2 - m\ + ie) {q 2 - m\ + ie) ' 1 J 

where 

e a = (-l) 1+ ^. (A6) 

The mass mx is now defined as m 2 x = | {n^jj + m u) ■ The flavor neutral propagator can be 
conveniently rewritten as 



where 



^VaVb ^-a&abPa + e a e b H ab (P a , P b , P x ) , (A7) 



p _ p p 

1 a n n . 5 1 b 9 ry . ) 1 X 9 2 I • ' 

q z — mf ja + ie q z — m* b + ie q z — m x + ie 



U ab (A,B,C) 



(ml- 



jj 



[m\ - m)j) (m 2 x - m 2 



m 



•la 



m 



>ib 



)( 



m 



lb 



m 



x 



B 



m 



x 



m 2 J (m 



x 



m. 



2 ) 

1/6 / 



c 



(A8) 



To include the baryons with one heavy quark into the theory, we use the same interpolat- 



ing fields, Eqs. (fTij) and (fT7 |l . The s e 



The baryon tensor of QCD is con- 
2. In our normalization 



baryons are then described by the field %j which 
forms a 17-dimensional representation of S77(4|2) 
tained as %j = T^, when the indices are restricted to the range 
convention, we have 

J_ ( Aq'' 
V2\~Aq 



T ■ 



(A9) 



The 



Se — 1 baryons are described by which makes up a 19-dimensional representation 
of SU{4\2). The baryon tensor of QCD is embedded as Sf, = S^, when the indices are 
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TABLE VI: The coefficients and A^, in S77(4|2) PQ^PT. Coefficients are listed for the Aq, 
and for are grouped into contributions from loop mesons with mass m^, while for A^i are 
grouped into contributions from pairs of quark-basis r\ q mesons. 
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restricted to the range 1 — 2. Here the QCD flavor tensor S^- is given in terms of Bij and 
B*j as in Eq. (HI, but with 

/ S+ 1 -^£° \ 

^ = 1 2,o > ( A1 °) 

and similarly for B*j*. The superscript on these states labels the 3-projection of isospin. The 
remaining states relevant to our calculation have been classified in [27] . 

The free Lagrangian for the T and fields is the same as Eq. (121)) and the Lagrangian 
that describes the interactions of these fields with the pseudo-Goldstone modes is given 
by Eq. (J2EJ). The LECs appearing in the Lagrangian, the Xj, #2, and g^, are of course 
numerically different than those in S77(6|3). The constants Ai, A2, gi and #3 have the same 
numerical values as those used in SU(2) heavy hadron xPT. The remaining constants A!j P< ^ 
and X^f^ are different in the partially quenched theory and the LEC A4 of SU (2) xPT can 
be found from matching, see Appendix [Bl 



2. Baryon masses 



For the mass of the Aq baryon in S77(4|2) PQxPT, the tree-level operators give the 
leading-order contributions, namely 



M 



(i) 



^Af Q) tr(mA + Ai PQ) strK). 



(All) 



The loop diagrams depicted in Fig. ^ make contributions to the Aq mass at NLO in the 
chiral expansion. Explicitly we find 



M 



(3/2) 



gi 



4tt 2 / 2 



2j A^K, A, 11) + 2J A^T{m^ mp, A, //) 



(A12) 



The non-analytic function T(m, A,/i) is given in Eq. (j3*T|) . while the hairpin loop function 
J r (m<f > ,m ( f ) r, A, ft) appears in Eq. (j32J). The coefficients and A^y in the sums over loop 
mesons and pairs of flavor-neutral mesons are listed for the Aq in Table IVII 

For the masses of the sg = 1 baryons in S77(4|2) PQ^PT, we have the tree level contri- 
butions 

M$ = X 1 m B + A 2 str(m ff ), (A13) 

where the coefficients itlb for non-degenerate quarks are listed for individual 5-states in 
Table EH 
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TABLE VII: The tree-level coefficients in xPT and PQxPT. The SU(2) and SU(A\2) coefficients 
mg are listed for the S£ = 1 baryon states B. 



B 


m B 




m u 




\{m u + m d ) 







TABLE VIII: The coefficients A$ and A^, in SU{A\2) PQ^PT. Coefficients are listed for the 
S£ = 1 baryons B, and for A B are grouped into contributions from loop mesons with mass m^, 
while for A^,, are grouped into contributions from pairs of quark-basis Vq mesons. 
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The loop diagrams that contribute to the mass are depicted in Fig. |21 and these include 
hairpin interactions. At NLO, we find 



M 



(3/2) 



gi 



127T 2 / 2 



9l 



12nf 



<f> tjxj,' 



(A14) 



The non-analytic functions T are given in Eqs. flUD and flH2J). The coefficients and Ajfc 
in the sums over loop mesons and pairs of flavor-neutral mesons are listed for the B baryons 
in Table IVILT1 and similarly for B? and B?,, in Table HX1 



TABLE IX: The coefficients B£ and Bfy in SU{A\2) PQxPT. Coefficients are listed for the s £ = 1 
baryons B, and for B B are grouped into contributions from loop mesons with mass while for 
B?j, are grouped into contributions from pairs of quark-basis r\ q mesons. 
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TABLE X: The coefficients in SU(3) xPT. Coefficients are listed for the se = baryon states 
T. 
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APPENDIX B: HEAVY HADRON MASSES IN SU(3) AND SU{2) 



For completeness, we include calculations of the baryon masses in ordinary x?T. These 
were calculated to one- loop order in [33] - To compare with the main text, we consider the 
case of SU(3) flavor in the isospin limit and to compare with the results of Appendix El we 
consider non-degenerate 577(2). For these calculations, we retain the tensors Tjj and in 
%PT. For the case of SU(3) flavor the tensors are given in Eqs. (fTfij) and (}2Tj|) . respectively; 
while for SU(2) flavor, they are given by Eqs. (jA9|) and ([AlOjh respectively. 

The free Lagrangian for the and fields in SU(3) and 577(2) both have the form 

£ = -i (g»v ■ DS„) + A (S^S,) + Ai (S^M+S^) + A 2 (S^S,) tiM + 

+i (Tv ■ DT) + A 3 (TM+T) + A 4 (TT) trM + . (Bl) 

The only exception is that for the SU(2) flavor group, the term with A3 is not independent; 
hence, one can choose A3 = 0. When matching the S77(4|2) Lagrangian onto the SU(2) 

Lagrangian, one finds A4 = ^X^f + X^^. 

The interaction Lagrangian for T^- and has the same form for both SU(3) and SU(2), 
namely 

£ = 192 (S^APS*) W + V2g 3 [(TAPS,) + (SM^T)] . (B2) 

The factor v^2 is chosen so that g 3 agrees with the literature, e.g. [43]. The numerical values 
of the Aj, #2, and #3 are of course different for the SU(3) and SU(2) theories. 



1. SU{3) 

The masses of the si = baryons in three flavor xPT are given at LO by 

jl4 1} = A 3 m T + A 4 tr(m g ) , (B3) 

where the coefficients rrij- are listed for individual T-states in Table |I] At NLO one has 
contributions from the loop diagrams in Fig. Q with the exception that there is no hairpin 
interaction in %PT. The expression for the mass is thus 

M^ = -4j- 2 J2AlH^A^), (B4) 

71 1 <f> 

where the function T appears in Eq. (}3*Tj) and the coefficients A^ are listed in Table 1X1 
Next the masses of the se = 1 baryons in three flavor x?T are given at LO by 

il4 1} = X 1 m B + A 2 tr(m 9 ), (B5) 
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TABLE XI: The coefficients in SU(3) xPT. Coefficients are listed for the S£ = 1 baryon states 
B. 
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TABLE XII: The coefficients B^ in SU(3) xPT. Coefficients are listed for the sg = 1 baryon states 
5. 
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where the coefficients for quarks are listed for individual S-states in Table ITTT1 The loop 
diagrams that contribute to the mass are depicted in Fig. El excluding the diagrams with 
hairpin interactions. The NLO expression for the mass is 



M 



(3/2) 



9l 



\2ix 2 p 



9l 



Unf- 



(B6) 



where the function T appears in Eq. (pTTj) . and the coefficients are listed in Table IXTl 
while B? appear in Table 1X111 



2. 517(2) 

The mass of the Aq baryon in two flavor x?T is given at LO by 

M { l ] Q = A 4 tr(m 9 ). (B7) 

The loop diagram without hairpin interaction shown in Fig. ^ makes the NLO contribution 
which reads 

M^ = ^^,A,^, (B8) 

with the function T defined in Eq. (KUJ) . 

Finally the masses of the Sg — 1 baryons in two flavor x?T are given at LO by 

Mjp = Ai m B + A 2 tr(m g ), (B9) 

where the coefficients itlb for non-degenerate quarks are listed for individual 5-states in 
Table IVIII The NLO expression for these masses arises from the loop diagrams of Fig. |21 
excluding the hairpins. We find 

Mi B /2) = -^Sp^m,, -A, + (BIO) 
19 
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